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Abstract 

We study the stability and robustness of large-scale vehicular formations, in which each vehicle is modeled as a 
double-integrator. Two types of information graphs are considered: directed trees and undirected graphs. We prove 
stability of the formation with arbitrary number of vehicles for linear as well as a class of nonlinear controllers. 
In the case of linear control, we provide quantitative scaling laws of the stability margin and sensitivity to external 
disturbances (Hoo norm) with respect to the number of vehicles N in the formation. It is shown that the formation 
with directed tree graph achieves size-independent stability margin but suffers from high algebraic growth of initial 
errors. The stability margin in case of the undirected graph decays to as at least 0(1/N). In addition, we show 
that the sensitivity to external disturbances in directed tree graphs is geometric in k where k < N is the number 
of generations of the directed tree, while that of the undirected graph is only quadratic in N. In particular, for 1-D 
vehicular platoons, we obtain precise formulae for the Hoo norm of the transfer function from the disturbances to the 
position errors. It is shown that the H x norm scales as 0(a N )(a > 1) for predecessor-following architecture, but 
only as 0(N 3 ) for symmetric bidirectional architectures. For a class of nonlinear controllers, numerical simulations 
show that the transient response due to initial errors and sensitivity to external disturbances are improved considerably 
for the formation with directed tree graphs. However, by using the nonlinear controller considered, little improvement 
can be made for that with undirected graphs. 

Index Terms 

Multi-agent systems, Stability margin, Convergence rate, Hoo norm, String instability, Nonlinear control, Dis- 
tributed control. 

I. Introduction 

Distributed control of multi-agent systems has spurred an extensive interest in the control community because 
of its wide range of applications such as automated highway system [1], collective behavior of bird flocks and 
animal swarms [2], and coordination of aerial, ground, and autonomous vehicles for surveillance, energy savings, 
mine sweeping etc. [3], [4], [5]. A typical issue in distributed control is that as the number of agents in the system 
increases, the performance of the closed-loop degrades progressively. A classical problem in which this progressive 
loss of performance was observed is the coordination of a platoon of vehicles moving in a straight line, see [1], 
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[6], [7], [8], [9] and references therein. The goal in the platoon control problem is to maintain a constant separation 
between pairs of nearby vehicles using a control law that is distributed, meaning that each vehicle can only use 
measurements of relative position and/or velocity with respect to its (at most two) neighbors that can be obtained 
with on-board sensors. 

Two decentralized control architectures that are commonly examined in the literature are predecessor-following 
and bidirectional architectures. In the predecessor-following architecture, the control action on each vehicle only 
depends on the relative information from its immediate predecessor, i.e. the vehicle in front of it. In the bidirectional 
architecture, the control depends on the relative information from both its immediate predecessor and follower. 
Within the bidirectional architecture, the most commonly analyzed case is the symmetric bidirectional architecture, 
in which the control at a vehicle depends on the information from both of its neighbors equally. 

It has been established that the predecessor-following architecture suffers from high sensitivity to disturbances 
with linear control, see for instance [10], [11] and references therein. High sensitivity to external disturbance is 
typically referred to as string instability [6], [12] or slinky-type effect [8], [9]. Seiler et. al. showed that with 
linear control, string instability with the predecessor-following architecture is independent of the design of the 
controller, but a fundamental artifact of the architecture [7]. String instability can be ameliorated by non-identical 
linear controllers but at the expense of the control gains increasing without bound as the number of the vehicles 
increases [9], [13]. It was later shown in [7], [14], [15] that the symmetric bidirectional architecture, also suffers 
from poor sensitivity to external disturbances. 

The degeneration of performance occurs not only in 1-D vehicular platoons, but also in vehicular formations 
moving in higher dimensional space. Mesh stability, which is a generalization of string stability is proposed in [16] 
to study look-ahead interconnected systems. In [17], [18], it is shown that if the information graph used is connected 
and undirected, the maximum error can be made independent of the size of the formation only if there is at least one 
vehicle communicating with at least O(N) other vehicles in the formation. However, it was later shown in [18], if 
there are at least two vehicles communicating to O(N) other vehicles in the formation, the system will be unstable 
if its size is beyond a critical value. Moreover, if the information graph has bounded degrees, the maximum error 
due to disturbances will be at least O(N). 

Although a rich literature exists on sensitivity to disturbances for predecessor-following and symmetric bidirec- 
tional architectures with linear control, to the best of our knowledge, a precise quantitative measure and comparison 
of sensitivity to disturbances of these two architectures is lacking. Moreover, most of the work on formation control 
of platoons has been limited to linear control laws, while little is known about nonlinear control. Nonlinear terms 
in the closed loop dynamics may arise from either purposefully designed nonlinear control laws (if beneficial) 
or unavoidable non-linearities in the vehicle dynamics, such as actuator saturations. Both of these cases can be 
analyzed by considering linear plant dynamics and nonlinear controllers. 

In this paper we examine the stability and robustness (sensitivity to external disturbances) of vehicular formations 
with linear as well as a class of nonlinear controllers, for two types of information graphs: directed trees and 
undirected graphs. The 1-D vehicular platoons with predecessor-following and symmetric bidirectional architectures 
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are special cases of the formations we considered. Each vehicle in the formation is modeled as a fully actuated point 
mass (double-integrator). This is a commonly used model for vehicle dynamics in studying vehicular formations, 
which results from feedback linearization of actual non-linear vehicle dynamics [9], [19]. A few authors have used 
first order kinematic models by ignoring vehicle inertia. However, it is pointed out in [18] that kinematic models 
fail to reproduce phenomena such as high sensitivity to disturbance that are exhibited by large platoons with kinetic 
models of vehicles. We prove stability of the closed loop with an arbitrary number of vehicles for a class of non- 
linear controllers, where the control gain functions satisfy certain sector conditions. The difference between the 
transient responses with different type of graphs in case of linear control is explained by the expressions we derive 
for the least stable eigenvalue of the closed-loop state matrix and its multiplicity. In particular, we show that the 
formation with directed tree graph has a larger stability margin compared to that with undirected graph: constant 
vs. 0(1/ N). However, it suffers from algebraic growth of initial conditions due to the high multiplicity of the least 
stable eigenvalue. For the non-linear control, we study the transient performance through numerical simulations, 
which indicate that the transient response of the formation with directed tree graph can be improved significantly 
by using a saturation-type nonlinearity in the control. 

Next, we examine the closed-loop's performance in terms of the sensitivity to external disturbances. Specifically, 
we examine the amplification factor, defined as the L-2 gain from the disturbances acting on all the vehicles to 
their position tracking errors. In case of linear controllers, we prove that the amplification factor, which becomes 
a Hoo norm, grows as 0(a k ) {a > 1) for directed tree graphs where k is the number of generations of the 
directed tree, but only as 0(N 3 ) for the that with undirected graphs. For 1-D vehicular platoons, k = N in 
the predecessor-following architecture. This establishes a precise comparison between the symmetric bidirectional 
and predecessor-following architectures in terms of robustness to disturbances. In addition, we show that the peak 
frequency for the predecessor-following architecture is 0(1) while for the symmetric bidirectional architecture it 
is 0(1/N). Establishing such scaling laws in case of non-linear systems has so far been proved intractable. We 
therefore study the response in the non-linear case through extensive numerical simulations, with both sinusoidal 
and random disturbances as inputs, and estimate performance metrics from simulation data. We observe from these 
studies that, within the predecessor-following architecture, a nonlinear controller with a saturation-type nonlinearity 
performs better than the corresponding linear one. In the symmetric bidirectional architecture, the difference between 
the linear and nonlinear controller's performance is not that significant, though these conclusions are valid only for 
the specific controllers we investigated. 

In terms of the stability analysis with non-linear controllers, our work closely parallels that of [20]. However, the 
results of [20] only considers absolute velocity feedback while we consider relative velocity feedback. Furthermore, 
the assumption of symmetry made in [20] precludes the predecessor-following architecture from their formulation. 
In terms of sensitivity to external disturbances with linear control, our work is related to [18], [17], [21] and [22]. 
The discussion on [18], [17] have been stated above. In [21], it was shown that for a 1-D vehicular platoon, a 
disturbance sensitivity metric grows linearly in N for the symmetric bidirectional case. Veerman [21] also examined 
an asymmetric bidirectional architecture, though not the predecessor-following one, and showed that amplification 
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of errors grow exponentially in N in that case. Bamieh et al. [22] obtained scaling laws of certain H2 norms from 
disturbance to outputs that quantify a number of performance measures such as local error, long-range deviation 
etc. In this paper we consider the Hoo norm metric. 

The rest of this paper is organized as follows. Section II presents the graph theory and the problem statement. The 
stability analysis appear in Section III. The results on sensitivity to external disturbances are stated in Section IV. 
The paper ends with a conclusion in Section V. 

II. Graph theory and Problem statement 

A. Graph Theory 

In this paper, we use an information graph to model the interaction topology between controlled vehicles. An 
information graph is a graph G = (V, E) with vertex set V and edge set E. The set of edges E C V x V 
specify the information flow used to compute their local control actions. We consider two types of information 
graph: directed graph and undirected graph. In a directed graph, we use to represent a directed edge from 
i to j, in which i is called the parent of j and j is called the child of i. The set of neighbors of i is defined as 
Mi := {j £ V : (j, i) £ E}. In particular, we restrict ourselves to a special class of directed graphs, directed trees. 
A directed tree is a directed graph which would be a tree if the directions on the edges were ignored. We say a 
directed tree has a spanning tree if every node, except one special node without any parent, which is called the root, 
has exactly one parent, and the root can be connected to any other nodes through paths. For a directed tree graph 
who has a spanning tree, we call the nodes whose parent is the root (reference vehicle) as the first generation, the 
nodes whose parent is in the first generation as the second generation and so forth. Without loss of generality, we 
assume there are rik nodes in the fc-th generation, and the nodes in it are labeled as {rik-i + 1, • • ■ , rik-i + Hfc} 
by identifying iiq = 0. 

In an undirected graph, the edges have no direction, i.e. i £ Afj if and only if j £ M- We say an undirected 
graph has a spanning tree if there is a tree composed of all the vertices and some (or perhaps all) of the edges. For 
an undirected graph, it has a spanning tree if and only if it's connected. In both directed and undirected graphs, 
the graph Laplacian L is defined by = — 1 if e 3 -j £ E and La = — J^ fc L^. The grounded (Dirichlet) graph 
Laplacian L g of the information graph is obtained by removing from L the row and column corresponding to the 
grounded node (reference vehicle). 

B. Problem Statement 

We consider the formation control of TV + 1 homogeneous vehicles which are moving in Euclidean space. For 
ease of exposition, we only consider one dimension of the translation motion. The analysis is also applicable to 
all three dimensions, as long as the dynamics of each of the coordinates of a vehicle's position are decoupled and 
each coordinate can be independently controlled. Under this fully actuated assumption, the closed loop dynamics 
for each coordinate of the position can be independently studied; see [23], [18] for examples. The position of the 
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Fig. 1. Desired geometry and its information graphs of a vehicular platoon with N vehicles and 1 "reference vehicle", which are moving in 
ID Euclidean space. The filled vehicle in the front of the platoon represents the reference vehicle, it is denoted by "0". 



i-th vehicle is denoted by pi and each vehicle is modeled as a double integrator: 

pi=Ui+Wi, ie {1,2,- ■■ ,N}, (1) 

where Ui is the control input, and u>i is the external disturbance. This is a commonly used model for vehicle 
dynamics in studying vehicular formations, which results from feedback linearization of actual non-linear vehicle 
dynamics [9], [19]. 

The control objective is that vehicle maintain a rigid formation geometry while following a desired trajectory. 
The desired geometry of the formation is specified by the desired gaps A^j) for i, j S {1, • ■ ■ , N}, where A^^j) 
is the desired value of pi(t) — pj(t). The desired inter-vehicular gaps A^^Vs are positive constants and they have 
to be specified in a mutually consistent fashion, i.e. A-Uk) = + ^(j,k) f° r every triple (i, j, k). The desired 

trajectory of the formation is provided in terms of a reference vehicle. Since we are interested in translational 
maneuvers of the formation, we assume there is only one reference vehicle, which is denoted by index 0. The 
trajectory of the reference vehicle is denoted by Po(t) and its dynamics is independent of the other vehicles. The 
desired trajectory of the i-th vehicle, p*(t), is given by 

pJ(t)=pS(t)-A (0l< ). (2) 

In this paper, we consider the following decentralized control laws, where the control action on each vehicle 
only depends on the relative position and velocity information from its neighbors in the information graph: 

Ui = - ~ Pi + A (J.i)) + 9(Pi - Pji) ■ (3) 

In the above controller, we assume the possibly nonlinear functions /, g : M — > M are odd functions, which are 
smooth enough to guarantee the existence of solution of the coupled ODEs. Note that the information needed to 
compute the control action can be easily accessed by on-board sensors such as radar and doppler sensors. 

Moreover, in this paper, we are particularly interested in a special class of formations, i.e. vehicular platoons. 
Figure 1 shows a pictorial representation of a vehicular platoon and its directed and undirected information graphs 
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we considered in this paper. Corresponding to each type of the information graph (directed tree and undirected 
graph), we consider the following two decentralized control architectures: 

1) Predecessor-following architecture. The control action at the i-th vehicle depends on the relative measurements 
from its immediate front neighbor, which results in the following control law: 

Ui = - f(Pi - Pi-i + - 9(pi - Pi-i)- (4) 

2) Symmetric bidirectional architecture. The control action at the i-th vehicle depends equally on the relative 
measurements from its immediate front and back neighbors, which results in the following control law: 

Ui = - f{Pi - Pi-i + - g(pi - Pi-i) 

- fip% - Pi+i - A ( . M+ i)) - g(pi - Pi+i), 

UN = - f(pN - PN-l + A(JV-1,JV)) - g(PN - PN-l), (5) 

where i e {1,2,- •• ,N-1}. 
For a vehicular formation which is desired to translate at a constant speed, we assume the reference trajectory is 
to be a constant velocity type, i.e., Po(t) = Vat + cq for some constants Vq, Co- To facilitate analysis, we define the 
following position tracking error: 

p 2 := p t -p*, (6) 

where p* is given in (2). The closed-loop dynamics of the vehicular formation can now be expressed as the following 
coupled-ODE model 

Pi = ~ 5^ (/(Pt ~Pj) + 9 (Pi ~ Pj)) + w ti (7) 

where i € {1,2, ••• ,N}. Note that po(t) = po(t) = since the reference vehicle perfectly tracks its desired 
trajectory. The system can be expressed in the state space form: 

x = f (x, w), (8) 

where the state and disturbance vectors are defined as x :— \px,pi, ■ ■ ■ ,Pn,Pn] T and w := [w\, • ■ ■ , wm] T ■ The 
special case f(z) = k$z and g(z) = b$z (where ko > 0,&o > 0) in the above coupled-ODEs correspond to the 
case of linear control. In that case, the closed-loop can be represented as: 

x = Ax + Bw, (9) 

where A is the state matrix that depends on ko, bo and B is the input matrix with appropriate dimension. 

III. Stability Analysis 

In this section, we present the stability analysis of the undisturbed vehicular formations with both linear and 
nonlinear controllers. For the linear case, we also derive formulae showing the relationship between the stability 
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margin and the smallest and largest eigenvalues of its grounded graph Laplacian. The stability margin is defined as 
the absolute value of the real part of the least stable eigenvalue of the state matrix A in (9). The stability margin 
quantifies the system's convergence rate with respect to initial errors. For the case of non-linear controller, we 
provide sufficient conditions for asymptotic stability. Since convergence rates for non-linear systems are difficult to 
obtain analytically, we perform numerical simulations to study the convergence rate and transient performance with 
non-linear controllers and compare with that with linear controllers. In particular, we take 1-D vehicular platoon as 
examples. All simulations for studying transient performance correspond to the following scenario: we perturb the 
initial position of the first vehicle in the platoon from its desired value and observe the position tracking error of 
the last vehicle pN(t). In addition, for the convenience of comparison, we define the following energy measure of 
transient performance: 

i r T 

E := lim -j / p 2 N (t)+p%(t) dt, 

where xq is the initial error of the first vehicle, i.e., pi(0) = xq. We assume the above limit exits , i.e. the last 
vehicle has finite £2 energy. In numerical simulations, we use the following estimate of E, 



E 



1 



b o Jo 



PtfW+PArW dt, 



(10) 



where T is sufficiently large such that all the errors die out. We study through numerical simulations how E scales 
with the number of vehicle N and the initial error xq. 



A. Stability analysis with linear control 

With linear controller f(z) = koz, g(z) = boz (kg > 0, 60 > 0) in (3), it's straightforward to see that the state 
matrix A can be expressed in the following form, 

A = I N <g> Ai +L g <g> A 2 , (11) 

where ® denotes the Kronecker product, Tjv is the N x N identity matrix and A\ , A2 are given as below 

(12) 

The matrix L g is the grounded (Dirichlet) graph Laplacian of the vehicular formation with corresponding information 
graph. For example, the ground graph Laplacian for the 1-D vehicular platoon with the directed tree (predecessor- 
following) and undirected (symmetric bidirectional) information graphs are given as below: 
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(13) 
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We now present a formula for the stability margin of the formation in terms of the eigenvalues of its grounded 
graph Laplacian. 

Theorem 1: Consider a vehicular formation with linear controller, whose state matrix A is given in (11). If the 
information graph has a spanning tree, then the system is globally exponential stable and its stability margin is 
given by 



bgXi 
2 



2kg 



I I 2 fco- 



if X N < 4k /b 2 Q , 
if Ai > Ak Q /bl 
= |, otherwise. 



(14) 



where < Ai < • • • < A^v are the eigenvalues of the grounded graph Laplacian L g . Furthermore, 

1) With directed tree graph, the stability margin is S 1 -^ = 0(1) and the least stable eigenvalue of the closed-loop 
state matrix A is s m i n = bn+ V b o _ Ijfg^ m j s eigenvalue occurs with multiplicity N. 

2) With undirected graph with bounded degree, when N is large, the stability margin satisfies < £ for 
some positive constant c*. □ 

Proof of Theorem 1. Our proof follows a similar line of attack as [21, Theorem 4.2]. From Schur's triangularization 
theorem, every square matrix is unitarily similar to an upper-triangular matrix. Therefore, there exists an unitary 
matrix U such that 

U~ l L g U = L u , 

where L u is an upper-triangular matrix, whose diagonal entries are the eigenvalues of L. We now do a similarity 
transformation on matrix A. 

A ^{U" 1 ® I 2 )A(U®I 2 ) = (U^ 1 <g) I 2 )(In ® Ai + L g ® A 2 ){U ® 7 2 ) = In ® A 1 +L U ^A 2 . 

It is a block upper-triangular matrix, and the block on each diagonal is A\ + XeA 2 , where \( £ a(L g ), where cr(-) 
denotes the spectrum. Since similarity preserves eigenvalues, and the eigenvalues of a block upper-triangular matrix 
are the union of eigenvalues of each block on the diagonal, we have 



a(A) = a(A)= (J {a(A 1 + X e A 2 )} = (J {a 

X e £a(L g ) \ e £a(L g ) 

It follows now that the eigenvalues of A are the roots s of 

s 2 + A^feo s + Affco = 0. 





-fc Af 



1 

-b X e 



(15) 



(16) 



We now prove that the eigenvalues A/s are real and positive. For a directed tree graph who has a spanning tree, 
by the labeling method given in Section II-A, it's straightforward to see that the ground graph Laplacian L g d ^ is 
a lower triangular matrix whose diagonal entries are all l's. Thus all its eigenvalues are real and positive. For 

(u) 

a connected, undirected graph, it follows from [24, Lemma 1] that the ground graph Laplacian L g is positive 
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definite. In addition, it's symmetric, so its eigenvalues are real and positive. Since fco > 0, 60 > 0, Xi > 0, we obtain 
from (16) that all the eigenvalues of A has negative real part. Therefore, the system is globally exponential stable. 
For each I € {1, 2, ■ ■ ■ , N}, the two roots of the characteristic equations (16) are denoted by sf, 



S ^~— ± 2 • 

The one that is closer to the imaginary axis is denoted by sf, and is called the less stable eigenvalue between the 
two. The least stable eigenvalue is the one closet to the imaginary axis among them, it is denoted by s m ; n . 
Depending the discriminant in (17), there are three cases to analyze: 

1) If An < 4feo/6g, then the discriminant in (17) for each I is non-positive, then recall that the stability margin 
S is defined as the absolute value of the real part of the least stable eigenvalue, which yields 

S = |i?e(s m in)| = -y^-- 

2) If Ai > 4fco/6p, then the discriminant in (17) for each I is non-negative, the less stable eigenvalue can be 
written as 

+ _ X e b Q - y/ (X(b ) 2 - 4X e k 2ko 

The least stable eigenvalue is achieved by setting A^ = Ajv, then have the convergence rate 

2k Q 



S = |i?e(s n 



60 + V b o ~ ik o/X N 

3) Otherwise, if the discriminant in (17) is indeterministic, i.e. it's negative for small I and positive for large I, 
then the stability margin is given by taking the minimum of the above two cases. 
In the case of directed tree graph, we know that the eigenvalues of the grounded graph Laplacian L g with 
directed tree graph are Ai = • • • = Xn = 1. From (17), we see that the eigenvalues of the state matrix A are 



(— 60 ± ^/£>g — 4fco)/2 and they have multiplicity of N, They are independent of N, thus the stability margin 
is O(l). For an connected, undirected information graph whose maximum degree is a function q(N + 1), where 
N + 1 is the number of vehicles (including the reference vehicle) in the formation. It follows from Lemma 3.2 
of [17] that Ai < q(N + 1)/N. If q(N + 1) is bounded, it follows from the three cases discussed above that when 
N is large, S^ u ' < c* /N for some positive constant c*. ■ 

Remark 1: For the connected, undirected graph, the upper bound on Ai can be tightened if the information graph 
is known beforehand [17], [25]. In particular, for 1-D vehicular platoons, we obtain from Theorem 2 of [26] that 
the eigenvalues Xi of L g u ^ are given by X( = 2 — 2 cos( - §^q^p ) = 4 sin 2 ( ^n+i) )• They are distinct. For large 
N and small I, we use the following approximation A^ = ^ — — , it's not difficult to see that the discriminant 
in (17) is negative. Therefore, the least stable eigenvalue is determined by setting I = 1, it has multiplicity 1. The 
smallest eigenvalue Ai is now 0(l/N 2 ), we obtain the stability margin S (u *> = 0(1/N 2 ). □ 

Although stability guarantees that transients due to initial conditions decay to as t — > 00, the convergence rate 
and transient performance depend quite strongly on the architecture. For a stable LTI system x = Ax, the stability 
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margin is an appropriate measure of this convergence rate. It follows from Theorem 1 that the stability margin 
of the formation with linear control decays to with increasing N for undirected information graph, while it is 
independent of N for directed tree graph. This makes the formation with directed tree graph architecture have faster 
convergence rate than that with undirected graphs, especially for large N. However, the large algebraic multiplicity 
of the least stable eigenvalue of the formation with directed trees will cause large algebraic growth of the initial 
errors before they decay to 0. For example, in a directed tree graph with k generations, the transient is proportional 
to £ fc e -ReOmm)t p or a directed tree graph with bounded degree, k — > oo as N — > oo. So for large N, the algebraic 
growth is obvious. Corroboration through numerical simulations is provided in Section III-C. 

B. Stability analysis with non-linear control 

The next results are on the stability of the formation with non-linear controllers. In the statements of the theorems 
that follow we say that a scalar function / belongs to the sector [e,K] if ez 2 < zf(z) < Kz 2 ,V z G K, and it 
belongs to the sector (0, oo] if zf(z) > 0, Vz^O. 

Theorem 2: Consider a vehicular formation with directed tree information graph. If it has a spanning tree, and 
f(z),g(z) satisfy the sector conditions / G [e\,Ki\,g G [£2,^2], where < E\ < K\ < 00, < £2 < K 2 < 00, 
then the origin x = is globally asymptotically stable. □ 

Theorem 3: Consider a vehicular formation with undirected information graph. If it has a spanning tree, and 
f(z),g(z) satisfy the sector conditions / G (0,oo],g G (0, 00], then the origin x = is globally asymptotically 
stable. □ 

Remark 2: Note that the stabilities of the linear controllers are special cases of Theorem 2 and Theorem 3. 
However, the latter two theorems have no implications on the stability margin. Comparing the above two theorems, 
we notice that the requirement on the sector condition with directed information graph is stricter than that of 
undirected information graph. However, we should note that these sector conditions are only sufficient. □ 

The proof of Theorem 2 will use the following proposition, whose proof is given in the appendix. 

Proposition 1: Consider the second order autonomous system 2/1 = 2/2,2)2 = —f{yi — ui) — 5(2/2 — 112), where 
2/1, 2/2, Mi, U2 G K and the odd functions f,g : R — > K lie in the sectors / G [ei, Ki], g G [£2,-^2], where 
< £\ < K\ < oo,0 < £2 < K2 < 00. The origin of the unforced system (with u(t) = [ui(t), U2(t)] T = 0) is 
globally exponentially stable (GES) and the system is input-to-state stable (ISS) with u as the input. □ 

Proof of Theorem 2. We observe that in a directed tree information graph who has a spanning tree, the dynamics 
of each node (vehicle) only depend on the position and velocity of its parent, but not other vehicles. For example, 
let i be an arbitrary vehicle in the graph, we use to denote the position error of its parent, the closed-loop 
dynamics of i can be written as: 

§i = -m-^ +) )-g(Pi-ft ) )- 

According to Proposition 1, the origin of the above unforced dynamics (with u(t) = [p\^\p^] T = 0) is GES, 
and it's ISS with u(t) as its input. 
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We first consider the subsystem consisted of the root (reference vehicle) and the first generation, whose closed-loop 
dynamics can be written as below by using the fact [po,po] T = 0, 



Pi 



-f(pi) - g(pi), 



,(m) 



{Prn =-f(p ni )-g(p ni )- 

where n\ is the number of nodes in the first generation and x^ ni ' = \pi,Pi,--- ,Pn 1 ,Pn 1 ] T ■ We notice from 
Proposition 1 that the origin of each dynamics in the above equation is GES, thus the origin of the cascade system 
x( ni > = f ni (a;( ni )) is also GES. We next consider the subsystem consisted of the root and the first two generations. 
Its closed-loop dynamics can be written as: 

Pi = -/(pi) - g(pi), 



Pm = -/(Pm) -g(Pm), 

Pm+i = -/(Pm+i - pi|+i) - g(Pm+i - p£+i), 



_(fii+n 2 ) _ f (r r (ni+n 2 )\ 
1 ni+n2V )•> 



^Pn 1 +n 2 — f{Pn 1 +n 2 P ni +n 2 ) 9(Pn 1 +n 2 P ni +n 2 )- 

where x^ ni+n2 ^ = [pi,pi,-"" ,Pn 1 +n 2 ,Pn 1 +n 2 ] T ■ Notice that the parents of the nodes in the second genera- 
tion are a subset of the nodes in the first generation, i.e. {p^ +1 ,p^ +1 , ■ ■ ■ , pi^+ n2 , PnH-n 2 } ^ s a surj set of 
{Pi>Pii ' ' ' iPnuPni}' The above dynamics can be divided into two parts: 



Jm+n 2 ) — f ( T (n 1 +n 2 )\ 
b 1 n 1 +n 2 V- 1 ' )i 



T (n 2 ) _ f ( r (n 2 ) r (ni)\ 
L n 2 \ ^ ; ) 7 



(18) 



where = [p„ 1+1 ,p ni+1 , ■ ■ ■ ,p ni+n2 ,p ni+ri2 ] T . The unforced system £ ( ™ 2) = f n2 (x ( ™ 2) , 0) is given by 

( 

Pni+l = -/(Pni+l) ~ ff(Pni+l)j 

x(™ 2 > = f„ 2 ( a > 2 >,0), => J : 

k Pni+n a = -/(Pm+na) _ ffCPni+na)- 

Again, following the same argument as proving a^™ 1 ' = ^(a^™ 1 )) is GES, we have that the origin of the unforced 



system x 



(»2) 



f n2 (x^ n2 \0) is GES. In addition, the functions /, g are smooth enough functions, we obtain 



from [27, Lemma 4.6] that the forced system x*-™ 2 ^ = f„ 2 (x*-™ 2 ^, a;*-™ 1 -*) is ISS with re*-™ 1 -* as its input. We now 
invoke [27, Lemma 4.7], the origin of the cascade system x^™ 1+ ™ 2 ' = i ni +n 2 (x(" 1+n2 )) given in (18) is globally 
asymptotically stable. By a straightforward induction method, we prove that the origin x = of the whole formation 
x = f(x) (x = \pi,pi, ■ ■ ■ ,Pn,Pn] T ) is globally asymptotically stable. ■ 
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Proof of Theorem 3. For the vehicular formation with connected, undirected information graph, we consider the 
following Lyapunov function candidate, which is inspired by the one used in [20]: 

i=o jeNi j=o 
where x = \px,Px, ■ ■ • ,Pn,Pn] T - The derivative of V along the trajectory of (7) with W{ = is 

N N 

^=EE - - fc) + 2 E^ 

i=o jejv", i=o 

AT JV 

=E E f(Pi-pj)(Pi-p~j)+ 2 J2p~i(- E (f(Pi-Pj)+9(pi-Pj))) 

AT 

= ~EE^~ -Pj)< °- 

i=o jeJVi 

If V = 0, then we have pi = for all i, since g(z) satisfies zg(z) > 0, Va; 7^ and po = by definition. Asymptotic 
stability now follows from LaSalle's Invariance Principle. In addition, we have V(x) — ► 00 as ||x|| — > 00, i.e. V(x) 
is unbounded in each direction of the Euclidean coordinate of x. Therefore, the Lyapunov function V is radially 
unbounded, and we get global asymptotic stability. ■ 

C. Numerical comparison between linear and nonlinear controllers for transient decay 

Throughout this section, we consider the following specific linear and nonlinear controllers. The control gain 
functions f(z) and g(z) used in controllers (4) and (5) are given by 

Linear: f(z) = k z, g(z) = b z, 

Non-linear: f(z) = Bx tanh(7iz), g(z) = B2 tanh(72z), (19) 

where fc = 1,6q = 0.5, B x =5,71 = 0.2, B 2 = 5,72 = 0.1. The parameters have been chosen in such a way 
that the slopes of f(z) of g(z) near the origin are equal to fco an d bo, respectively. This is done to make the linear 
and non-linear cases comparable to some extent. Note that these f(z) and g(z) do not satisfy the global sector 
conditions assumed in Theorem 2, but only satisfy the sector conditions locally. However, the region in which they 
satisfy the sector condition can be made arbitrarily large by choosing sufficiently small ex and £2. 

In this section, we compare the convergence rate and transient performance between linear and nonlinear 
controllers through numerical simulations. Figure 2 (a) depicts the transients of the platoon with linear and nonlinear 
controllers for predecessor-following architecture. The algebraic growth for linear controller which is predicted by 
Theorem 1 is observed. We also see that the nonlinear controller has much smaller peak error than the linear 
controller. The transients of the symmetric bidirectional architecture are shown in Figure 2 (b). We see that (i) the 
performance of the non-linear case is similar to that of the linear controller, and (ii) the peak value of the error is 
much smaller compared to that in the predecessor-following architecture. 
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} 100 200 300 1,5 100 200 300 

time (t) time (t) 



(a) Predecessor-following (b) Symmetric bidirectional 

Fig. 2. Comparison of transients of the position tracking error of the last vehicle for a platoon of TV = 10 vehicles between linear and nonlinear 
controller. The initial condition of the first vehicle used is xo = 10. 



Moreover, we see from Figure 3 that the convergence rates of the linear and non-linear controllers in the 
predecessor-following architecture are similar. In addition, the error in the predecessor-following architecture is 
smaller than in the case of symmetric bidirectional architecture for large t, irrespective of whether the control is 
linear or non-linear. This is due to the fact that in the predecessor-following architecture, the stability margin is 
independent of the number of vehicles TV in the platoon, while it's very small when N is large for symmetric 
bidirectional architecture. 

Figure 4 and Figure 5 show the estimate of energy measure E for T = 10 4 seconds (defined in (10)) as a function 
of N and xo respectively. We see that (i) the energy in the predecessor-following architecture has a much worse 
scaling trend with N or xq than that in symmetric bidirectional architecture, no matter the controller is linear or 
nonlinear, (ii) nonlinear controller performs better than linear controller in the predecessor-following architecture, 
whereas it performs similarly or worse in the symmetric bidirectional architecture. 

In summary, the predecessor-following architecture has a faster convergence rate than that of symmetric bidi- 
rectional architecture, no matter the controller is linear or nonlinear. However, due to the high algebraic growth 
of initial errors, the predecessor-following architecture has a worse scaling trend of energy, comparing to the 
symmetric bidirectional architecture. Moreover, if only one architecture is considered, the nonlinear controller in 
general has a better (respectively, worse) transient performance than that of linear controller in the predecessor- 
following (respectively, symmetric bidirectional) architecture. 
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Fig. 3. Comparison of transients of the position tracking error of the last vehicle for a platoon of N = 10 vehicles between predecessor-following 
and symmetric bidirectional architectures. The initial condition of the first vehicle used is xo = 10. 



IV. Sensitivity to external disturbances 

In this section, we study the sensitivity of the vehicular formation to external disturbances. Specifically, we 
examine appropriate gains from disturbances acting on all vehicles w <G K to the position tracking errors p £ M. N . 
We define the amplification factor AF as the L2 gain from the vector of disturbances w(t) = [wi(t), ■ ■ ■ , Wjv(i)] 
to position tracking error vector p(t) = • • • ,PN(i)]' 

^-plinear or nonlnear g^p Ijj^ij (20) 

IHUa ' 

In this paper, the £2 norm is well-defined in the extended space H\ = {u|u T £ £2,V r G [0, 00)}, where 
u T {t) = (i) u(t), if < t < t; (m) 0, if t > t. See [27, Chapter 5]. With a little abuse of notation, we suppress 
the subscript and write £2 = £|- 

In the linear case the amplification factor is equivalent to the Hoo norm of the transfer function G(s) from w to 

P. 

AF Unear = max (J max (G(jtj)) = (7 max (G(j Wp )), 
w 

where we have assumed the maximum is achieved, and u p := argmax w u max (G(juj)), cr max denotes the maximum 
singular value. In the non-linear case, evaluating AF nonhnear is intractable, so we use following conservative 
estimate: 

- nonlinear P f^M 

AF := , ( } , (21) 

\\w\\c 2 (t) 
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Fig. 4. Comparison of E between linear and nonlinear controllers as a function of N. The measure E is estimated by numerically evaluating 
the integral in (10) for T = 10 4 s. The initial condition of the first vehicle used is xq = 10. 



where ||u(i)||£ a ( T ) = y Jq J2iuf(t)dt, r is large enough and w = [a\ sm(oj p t + 6\), ■ ■ ■ ,ajy sin(uj p t + 9^)]. 
The parameters a = [»!,••• , Ojv] and 6 = [0±, - ■■ ,9n] are those that yield the £ 2 g am m the linear case. The 
choice of these parameter will be given later. Note that the estimates for the non-linear case are lower bounds: 

-s nonlinear , ■ , - 

AF < AF n * near . In the numerical simulations, the quantity AF ar for linear controller are also 

~ linear 

estimated for large r to compare with analytical results, which will be denoted by AF , as defined in (21). 

We also examine the effect of random disturbances. Specifically, let w(t) in the closed-loop platoon dynamics (8) 
be a vector of white noise with autocorrelation matrix aol, where ao is a constant and / is the identity matrix with 
appropriate dimension. Similar to sinusoidal disturbances, we define the following metric 



j^linear or nonlnear hni (22) 

t— too <7q 

where E(.) denotes the expected value and we have assumed the above limits exist. Notice in the linear case, the 
above ratios are exactly the H2 norms of the appropriate transfer functions from the white noise disturbances to 
the position tracking errors. The steady-state covariance matrix of the state p(t) of the system (9) that is driven by 
a white noise process w(t) is given by solution P of the following Lyapunov equation [28, Chapter 4]: 

AP + PA T = -Q, 

where Q = E[Bww T B T ], and B is the input matrix. Since A is Hurwitz, it guarantees the limit in (22) exists [28]. 
The steady-state expectations E(p(t) T p(t)) can be obtained by summing the odd diagonal entries of P, which 
yields 



v '££i^(2i-1.2i-l) 

j^linear V v | ^ (23) 

CTo 
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Fig. 5. Comparison of E between linear and nonlinear controllers for a platoon of N = 10 vehicles with different initial conditions xn- The 
measure E is estimated by numerically evaluating the integral in (10) for T = 10 4 s. 



For the non-linear controllers as well as linear controllers, we use the following estimate of the ratio defined in (22), 



■^linear or nonlnear P(--Q) (24) 

CTO 

where T is sufficiently large. Monte-Carlo simulations are used to estimate the ratios. For example, to estimate the 
ratio Ji nonhn ^ ar f or the 1-D vehicular platoon with predecessor-following architecture, the noise-driven system is 
converted into a standard stochastic differential equation (SDE) form 

dpi = pidt, dpi = -f(pi - Pi-i)dt - g(pi - pi-i)dt + a dWi(i), (25) 

where W(t) = [Wi(t), ■ ■ ■ , Wjv(t)] is a standard TV-dimensional Wiener process. Sample paths of the states are 
computed by using Euler-Maruyama Method to numerically integrate the SDE (25) [29]. The ratio Ji nonlnear i s 
now estimated by performing appropriate averaging over a large number of simulations, after letting each simulation 
proceed sufficiently long to allow transients to die out. 

A. Sensitivity to disturbance with linear control 

As stated earlier, analytical results on the sensitivity to disturbance is possible only for the linear case. The first 
result is on the sensitivity of the vehicular formation with directed tree graph. 

Theorem 1: Consider a vehicular formation with directed tree information graph who has a spanning tree. Let 
k be the number of generations. With linear controller f(x) — k$x and g(x) = box in (3), the amplification factor 

AF ltnear sat j sfies 

,fe-l / a rplinear ^ fh{& ~ 1) 



fta < AF ar < 



a 
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where a = \T(Jur)\ > 1, f}\ = \S(ju)r)\ an d $1 — \S(jws)\> m which 

TV A & OS + fco s 1 



+ 6 s + k ' s 2 + 6 s + fco ' 

and u>t and ujs are the peak frequencies of T(s) and S(s) respectively. □ 
Notice that in the above theorem, if the directed tree has bounded degree d, then we have k — > oo as N — > oo. 
The amplification factor is unbounded for increasing N. 

Proof of Theorem 1. The proof is inspired by the proof of Lemma 1 in [7]. With linear controller f(x) = k^x 
and g(x) = box, taking Laplace transform of the coupled-ODE model (7) and assuming zero initial conditions, we 
obtain the transfer function from the disturbances w — [wi, . . . , wn] t to position errors p = [pi, . . . ,Pn] T 

G(s) = (s 2 I N + (b s + k )L^)- 1 , 

where Zjv is the N x N identity matrix and L g is the grounded graph Laplacian of the directed tree. It follows 
from 1 1 G| | max < ||G||2 of [30, Chapter 10] that the Hoo norm of the transfer function G has the following relation 

AF Unear = || Gj |^ > \\G(jU T )\\ 2 > \\ G{ju T ) ||max = max |Gy (juJ T )\ • (26) 

We now look at a particular entry of G(s). Since there is k generations in the directed tree, let (0,pi), (pi,P2), ■ ■ • , (jPk-i-,Pk) 
be a longest path starting from the root. Its closed-loop dynamics can be written as 

P Pl = -k (p Pi -Ppi-i) ~ b (Ppi -Ppi-i) +W Pi , *£{!,'■■ 

Taking Laplace transform of the above coupled-ODEs and assuming zero initial conditions, we have the transfer 
function from the disturbance w pi on the first vehicle to position tracking error p Pk of the last vehicle 

G ^( s ) = W^= ( (b 2°lt k T,\ k - S ( s ns) k -\ (27) 
W pi (s 2 + b s + k ) k 

where P Pk and W Pl are the Laplace transforms of p Pk andw Pl respectively. In addition, \G WlPk (ju)T)\ < niaxj^- \Gij(ju)r)\, 
we obtain from (26) a lower bound for AF lmear 

AF Unear >\G WlPk {jw T )\ = f3 l a k -\ 

where a = \T{jojT)\, Pi = \S(jujt)\- In particular, straightforward algebra shows that a > 1. 
In addition, for any matrix A, its spectral radius satisfies p(A) < \\A\\i, we have 



l^lkoc = sup \\G(jcj)\\ 2 = sup y/p{G*{ju)G{ju))) < sup y/\\G*{jw)G{ju))\\ 



< sup ^WWMGWWi = sup V||G0'w)||oo||G(j^)|| 



< sup V||G0w)||oo sup j\\G{jw)\\i. (28) 

UJ UJ 

As the particular entry G WlPk (s) given in (27), each entry of G(s) is either or S(s)T(sy~ 1 for some 1 < j < k. 
The absolute row sum of G(s) is J2j \Gij(J u )\> w h ere Gij(jui) is the transfer function from disturbance Wj to 
pi. Because of the special coupling structure of the directed tree, Gij(jui) — if j is not a predecessor of i and 
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Gij(jtu) = S(s)T(s) 1 1 if j is the last (I — l)-th predecessor of i, i.e. if j = i, then I = 1; if j is the parent of i, 
then / = 2 and so forth. Since there is k generation in the directed tree, the maximum absolute row sum is given 
by 

k 

IIGtfaOlU = max^ \G tj (ju>)\ = £ ^(jo.) | |r(io;) | ^ . 
j i=i 

In addition, we have su PtJ 1 | ^(o;)^- 1 1 < sup w |5(jw)| su Pu IT^'- 1 ] = where /3 2 = |S(jws)|. 

Thus we have 

sup||G0 W )|| oo =sup^|5(^)||T(^)|^ 1 <^/3 2 a i - 1 = /3 2 — i (29) 
w ;=i i=i 

Following the same argument as before, we obtain that 

k i 

sup || GO") Hi <fo^—- ■ (30) 

U OL — 1 



i n ea r 



Substituting (29) and (30) into (28), we obtain a upper bound for AF 

a k - 1 

||G|| ffoo </3 2 -. ■ 

a — 1 

Corollary 1: Consider an iV-vehicle platoon with predecessor-following architecture. With linear controller f(x) = 
k$x and g(x) = box in (4), the amplification factor AF hnear satisfies 

fta"" 1 < AF Hnear < ^ - X \ 

a — 1 

Furthermore, when N S> 1, 

Af* Ba,r wj3 lJ^o -» ~ , -■ (31) 

y or — 1 &o 

Moreover, a sufficient condition for a disturbance w = [wi, • • • , u^] = [oi sin(cji + 0\), • • • , sin(u>t + On)] to 

yield the amplification factor is a = [oi, • • ■ , a^v] = [oi, 0, • • ■ , 0], where ai is an arbitrary constant and u = lu p , 

9=[6 1 ,--- ,9 N ]=0. □ 

Proof of Corollary 1. The inequality of A.F lmear follows immediately by noting that in the predecessor-following 
architecture, the information graph is a line graph who has N generations. 

To prove the asymptotic results, the Hoo norm of a system can be interpreted in a sinusoidal, steady-state sense 
as follows (see [31]). For any frequency cj, any vector of amplitudes a = [oi, • • • , ajv] and any vector of phases 
9 = • • • , 9n], the input vector 

w = [wi, ■ ■ ■ , wn] = [a i sin(wt + 0\), • • • , a,N sin(wi + On)] 

yields the steady-state response e of the form 

e= [pi,-- - ,pjv] = [&i sin(wt + V>i), ■ ■ • , 6jv sin(o;t + Vat)]- 
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The i/oo norm of G(juj) can be defined as 

||G(j W )|| ffoo = sup M^. ( 32) 

uieR+ ,a,9eR N \\ W \\C 2 

Because of the special coupling structure of the predecessor-following architecture, the disturbance acting on a 
vehicle only effect its descendants. In addition, the supremum of the norm of the complementary sensitivity function 
sup u |T(s)| > 1, therefore to get the worst amplification factor, the disturbance should all concentrate on the 
first vehicle, that is the other vehicles are undisturbed (w; = for i € {2, • • • , N}). In addition, the steady 
state response correspond to a sinusoidal disturbance on the first vehicle w\ = a\ sin(o;i) are given by pi(t) ~ 
a\\T(ju:)' t ~ 1 S(juj)\ sm(u>t + ipi). It follows (32), the Hoc norm of the system is given by 

[ ai \S(joj)\sm(ojt + • ' ' , ailTO^)*- 1 ^)! sm(ut + ip N )]\\c 2 



\G{ju)\\ Hoo = sup- 



|ai sin(a;£)|| 



sup , 



jr \T(pu,y-is(j U )\* = sup J\s(jw)\2 11 



Since a = sup w \T(juj)\ > 1, in the limit of N — >• oo, the supremum of the above equation is determined by 
sup w \T(jui)\ 2N , thus we have 



\\G(M\\h x =\\S(Jut 




\l \//?~+2fc 3 b 2 — k 2 

Following straightforward algebra, we have lj p = ojt = - bo - - — -, which completes the proof. ■ 

The next theorem is the corresponding result for the vehicular formation with undirected information graph. 
Theorem 4: Consider a vehicular formation with connected, undirected information graph. With linear controller 
f(x) = kox and g(x) = box in (3), the amplification factor and its peak frequency satisfy 



AF 



-t—j — , otherwise. | 0, otherwise. 



where Ai is the smallest eigenvalue of the grounded graph Laplacian 

Moreover, a sufficient condition for a disturbance w = [wi, • • • , wn] = [oi sin(wt + 6*i), • • • , a^v sm(u>t + 9n)] 
to yield the amplification factor is a = [ai, ■ ■ ■ , ajy] = v\, where v\ is the eigenvector of Lg corresponding to 
the principal eigenvalue Ai and u> = oj p and 9 = [6\, ■ ■ ■ , 9^] =0. □ 

Remark 3: The above result is complimentary to Theorem 3.1 of [17] and Theorem 6 of [18]. In [17], the 
authors considered undirected information graph whose maximum degree is a function q(N + 1), where N + 1 is 
the number of vehicles (including the reference vehicle) in the formation. It follows from Lemma 3.2 of [17] that 
•^l < <l(N + l)/N. From Theorem 4, the amplification factor is 0((N/q(N + l)) 3 / 2 ), which is the same as the 
result of Theorem 3.1 in [17]. In [18], the authors considered the undirected information whose Ai is 0(1/ N), 
and derived an lower bound of norm by letting the peak frequency to be 0. In the context of Theorem 4, it 
yields AF hnear = which implies AF lmear is at least 0(N). This coincides with the result of Theorem 6 

in [18]. □ 
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Proof of Theorem 4. It follows from straightforward algebra that the transfer function from the disturbance w = 
[wi, . . . , wn] t to position error p = [pi, . . . ,Pn] T is given by 

G(s) = ( S 2 I N + (b oS + k Q )L^)-\ 

(u) (u) 

where L g ' is the grounded graph Laplacian of the vehicular formation. Since L g ' is real and symmetric, there 
exists an orthogonal matrix P such that Lg 11 ^ = PAP T , where A = diag{Ai, • • • , Xn}- The transfer function now 
becomes 

G(s) = (s 2 / + (b s + k^L)- 1 = P 

G N (s)_ 

where Ge(s) = s ^ + i, \] s +\ e k ■ ^ can ^ e shown using straightforward calculus that for each eigenvalue Xg, the 
maximum amplitude and its peak frequency of Ge(s) are 

\ yir , h L , . a , if A/ < 2fe /6§, 
A t := max|G £ (jo;)| = I A '' W 4fc °-*< 6 o (33) 

[j^, otherwise. 

[ aAwMM , ifXe <2k /bl 
u>e := argmax \Ge(jw)\ = < (34) 

I 0, otherwise. 

It follows from straightforward algebra that A\ > A 2 > ■ ■ ■ > A^. 

The Hoo norm of G(s) is now given by 



\\G\\ Hoo = sup||G(jw)||a = sup y/X ma x(G*{ju)G(ju>)) 

1 

= sup max ■ 



- t v /(-w 2 + A f fc ) 2 + ^% 2 
= max max \Gi (juj)\ = max Af = A\. 

i u I 

where A\ is given in (33). 

To prove the second part of the theorem, we rewrite the coupled-ODE model (7) as 

p + boLp + koLp = V\ sm(ojit). (35) 

By the method of eigenfuction expansion [32], we can write p(t) = v ehe(t), where u/s are the eigenvectors 

of Lg. Substituting it into Eq. (35), we obtain 



y](vihi(t) + b Lvihi(t) + k Lv£h e (t)) = v% sin(wit). 



Due to superposition property of linear system, the above equation can be split into N ordinary differential equations 
by using Lvg = Xgvi, 

hi(t) + bo\ehi(t) + k Xeh 1 (t) = sin^i), 

h t (t) + b X e h e (t) + k Q X e h e (t) = 0, t e {2, • • • , N}. 
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Following straightforward algebra, the steady-state response of each ht(t) is given by 

hi(t) = Ai sin(wii + -0i), 

hi(t) = 0, te{2,---,N}, 
where A\ is given in (33). Thus the steady state response of p is given by p = v\Ai sin(c<Jii + ipi), which yields 

\\P\\c 2 _ _ Aplinear ^ 



\Mc 2 

Corollary 2: Consider an TV-vehicle platoon with symmetric bidirectional architecture. With linear controller 
f(x) — kox and g(x) = box in (5), the amplification factor satisfies 

( }=— ) (2TV + l) 3 < AF Hnear < ( 1=) (2N + l) 3 . 

Furthermore, when JV S> 1, the amplification factor and its corresponding peak frequency are asymptotically 

(u) 

Proof of Corollary 2. Following Theorem 2 of [26], the eigenvalues of the grounded graph Laplacian L g (given 
in (13)) of a 1-D vehiculear platoon are given by 

Ai = 2- 2cos( — - — ) = 4sin 2 (— — -). 

y 2N+l' V 2(2TV+1) ; 

It follows from the fact ^x < since < x,V.t <S [0, \ ) that 

(2Jvtl) 2 - Xe ~ (2N+ir (36) 
For any fixed ko, bo, when TV is large, we have Ai < 2kojb\, which implies 

M = 2 > 1 > 1 (2TV + 1) 3 

A? /2 6ov/4*o - Mbl ~ Xl /2 b Vh ~ b oVk~o~ ^ 
where the last inequality is obtained from (36). In addition, using Ai < 2kojb\ and < ( 2JV + X ) ; we obtain an 
upper bound 

A, = 2 < 2 < ( 2jV + 1 ) 3 . 

A? /2 6 \/4fco- Ai^ " A? /a 6oV5fo ~ 46 \/2fco" ' 

To get the asynptotic formula, when TV is large, we use the approximation Ai « jjf?- Therefore, Ai < 2kojb\ is 

2 

true for large enough TV irrespective of the values of ko and bo. Substituting Ai ~ jj^j into (33) and (34), we 
obtain that 

AF u near=Ai ^ &N 3 U u= Wl - V&r 



/Mott 3 ' P 2TV 
which completes the proof. 



August 9, 2011 



DRAFT 



22 



10' 



&;io 5 



10 u 



& Nonlinear (simulation) 

* Linear (simulation) 

o Asymptotic formula of Cor. 1 

Lower bound (Corollary 1) 

Upper bound (Corollary 1 ) 




44 44""* 



10° 



^ 10 4 



10' 



if Nonlinear (simulation) 

^ Linear (simulation) 

O Asymptotic formula of Cor. 2 

Lower bound (Corollary 2) 

Upper bound (Corollary 2) 




10 

N 



30 



10 

JV 



30 



(a) Predecessor-following 



(b) Symmetric bidirectional 



Fig. 6. Comparison of amplification factor of the platoon with linear and nonlinear controller for both the predecessor-following and symmetric 
bidirectional architectures. 



B. Numerical comparison of sensitivity to disturbance between linear and non-linear controllers 

In this section, we present results of numerical computations in order to obtain the sensitivity to external 
disturbances of the 1-D vehicular platoon with nonlinear controllers. These are then compared against the predictions 
on similar performance measures for the linear controllers that have been analytically obtained in the previous 
section. The performance measures are also numerically computed for the linear controllers in order to verify the 
analytical predictions. The controllers used are the ones given by (19) in Section III-C. 

Figure 6 shows the amplification factor as a function of N. The following observations are obtained: 1) The 
lower and upper bounds and asymptotically formulae derived in Corollary 1 and Corollary 2 are quite accurate; 2) 
Comparing Figure 6 (a) and (b) we see that the symmetric bidirectional architecture has a smaller amplification 
factor than the predecessor-following case when the controller is linear. However, when nonlinear controller is 
applied, the symmetric bidirectional architecture has a worse scaling trend than that of the predecessor-following 
case; 3) In the predecessor-following architecture, the growth of the amplification factor with respect to N is much 
slower with the nonlinear controller than with the linear controller. In the symmetric bidirectional architecture, there 
is little difference between the two controllers for this sinusoidal disturbance. 

To examine the effect of random disturbances, we compute the estimate of R that is defined in (24) for T = 3000 
seconds, through Monte-Carlo simulations for both linear and non-linear cases. Figure 7 shows R vs. N for a fixed 
do, the strength of the white noise, while Figure 8 shows R vs. ctq for a fixed N. The following broad conclusions can 
be made from these simulations: 1) disturbance amplification is far worse in the predecessor following architecture 
than in the symmetric bidirectional one, and this conclusion holds whether one uses a linear or non-linear control; 
compare Figure 7 (b) vs. Figure 7 (a) and Figure 8 (b) vs. Figure 8 (a). This shows that architecture has a more 
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Fig. 7. Comparison of the ratios R, R as a function of the number of vehicles TV with white noise disturbances. The value of <tq used is 1. 



profound impact on performance than linearity or non-linearity, and given a choice, symmetric bidirectional is 
preferable. Though these conclusions are for a specific non-linear control design, simulations with other non-linear 
control laws that are not reported here due to lack of space shows the same trend. These laws include (a) f(z) = z 3 , 
g(z) = bz, and (b) f(z) = sat(z), g(z) = bz; 2) If symmetric bidirectional architecture is indeed used, both the 
linear and non-linear control laws have almost identical robustness to disturbances. The only exception is when the 
strength of the white noise <jq is large, in which case the non-linear control law performs poorly compared to the 
linear one; 3) If the predecessor architecture is to be used due to other constraints, the non-linear control law is 
preferable since it has better robustness to disturbance than its linear counterpart; see Figure 7 (a) and Figure 8 (a). 

V. Conclusion 

We studied the stability and robustness of large vehicular formations with both linear and a class of nonlinear 
decentralized controllers. Under certain sector assumptions on the non-nonlinear control gain functions, we proved 
the stability of the formations with both the directed tree and undirected graphs with spanning trees. In the case 
of linear control, we analyzed the transient response due to initial conditions in terms of stability margin and the 
multiplicity of the least stable eigenvalues. We showed that, on one hand, although the formation with directed 
tree graph achieves a size-independent stability margin, it suffers from high algebraic growth of initial conditions 
because of the high multiplicity of its least stable eigenvalue. On the other hand, the stability margin of the 
formation with undirected information graph decays to at least as 0(1/N) if it has bounded degree. Moreover, 
we derived the scaling laws of the amplification factor with respect to the size of the formation. In particular, for 1-D 
vehicular platoons, we show that the amplification factor scales as 0(a N ) (a > 1) for the predecessor-following 
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Fig. 8. Comparison of the ratios R, R of a 10-vehicle platoon as a function of the standard deviation tro of the white noises. 



architecture but only as 0(N 3 ) for the symmetric bidirectional architecture. We examined the transient performance 
and sensitivity to external (sinusoidal or random) disturbances with nonlinear controllers by extensive numerical 
simulations. These simulations show that in case of the directed tree graphs, a class of nonlinear controllers perform 
better compared to the linear control, both in terms of transient performance and sensitivity to external disturbances. 
In the case of undirected information graphs, numerical simulations indicate that the nonlinear controller yields equal 
or worse performance compared to the linear one. 

We have assumed that each vehicle has the same open-loop dynamics and uses the same control law. In the linear 
symmetric bidirectional case, it was shown in [33] that heterogeneity in vehicle masses and control gains has little 
effect on the stability margins. The non-linear stability analysis in this paper can be extended to the case when 
vehicles have different masses and employs a different controller in a straightforward manner. However, lack of 
symmetry has a large effect, as the difference between directed tree and undirected graphs established here shows. 
In the case of 1-D vehicular platoon with linear control, a general asymmetric architecture that lies in between 
symmetric bidirectional and predecessor-following architectures has been examined in [33] and [21]. It was shown 
in [33] that such asymmetry can lead to improvement in the stability margin. Numerical simulations in [34] indicate 
that if the asymmetry is only in the velocity feedback terms, sensitivity to external disturbances is improved as 
well. In contrast, it was shown in [21] that if equal asymmetry is applied to both position and velocity feedback 
terms, sensitivity to external disturbances worsens. Analysis of stability with non-linear asymmetric control and 
sensitivity to disturbance with asymmetric control (linear or non-linear) are open problems. 
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Appendix 

Proof of Proposition 1. First, we consider the unforced system with state y = [yi,y2] T , whose dynamics are 

Vi = 2/2, 

(37) 

m = -/(yi) -5(2/2)- 



Consider the following Lyapunov function candidate: 

•■in 



V(y) = ^y T P V + 1 J ' f(z)dz, 



(38) 



where P 



1 1 



and 7 > max {1, — + {1+K2) }, which ensures that P is positive definite. From the Rayleigh 



} 7^ 

Ritz Theorem [27], we have the following inequality A m i n (P)||j/|| 2 < y T Py < A max (P)||?/|| 2 , where A m j n (P) > 
0, A max (P) > are the minimum and maximum eigenvalues of P respectively. This shows that V(y) is radially 
unbounded, and in addition satisfies the following inequality 

<( p )n 112 , f Vl tl \ a ^ x ™™( p )u 1,2 , 2^ A max (P)+7^i ll „2 



V(y)<'-^^\\yr +1 J f(z)dz < ■-^f^Wyr + ^vl < m ~ v { '— \\v\F, (39) 

where the second inequality follows from the fact that the function f(z) belongs to the sector [ei, Ki\. The derivative 
of V along the trajectory of (37) is given by 

V = y T Py + 1 f{y 1 )y 2 

= -2/1/(2/1) - 72/23(2/2) + 2/2 + 2/12/2 - 2/19(2/2) 

< -e lV l - ( 7 e 2 - l)y 2 2 + (1 + K 2 )\ yi \\y 2 \, 

< -\{£ivl + (7^2 - 1)2/1) - 5^12/1 - 2(1 + ff 2 )Mllto| + (7£2 - 1)2/1)] 

< ™(ei»? + (7^2 - 1)2/1) < -\ min{ei, (7^2 - 1)}||2/|| 2 , (40) 

where the second last inequality following from 7 > maxjl,^- + ^j^-}, since that leads to the conclusion, 
upon a completion of squares, that Siyf — 2(1 + | £/i H2/2I + (7^2 — l)2/ 2 > 0. Since V is radially unbounded and 
satisfies (39), it follows from (40) that the origin y = of (37) is globally exponentially stable. Since the functions 
/, g are assumed to be smooth enough, the ISS property follows from the fact that a globally exponentially stable 
system with input u is ISS [27, Lemma 4.6]. ■ 
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